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Abstract. 

We present results from three-dimensional non-linear hydrodynamic simulations of a precession 
driven flow in cylindrical geometry. The simulations are motivated by a dynamo experiment 
currently under development at Helmholtz-Zentrum Dresden-Rossendorf (HZDR) in which the 
possibility of generating a magnetohydrodynamic dynamo will be investigated in a cylinder filled 
with liquid sodium and simultaneously rotating around two axes. 

In this study, we focus on the emergence of non-axisymmetric time-dependent flow structures 
in terms of inertial waves which - in cylindrical geometry - form so-called Kelvin modes. For a 
precession ratio (Poincare number) Po = Dp/Hc = 0.014 considered by us, the amplitude of the 
forced Kelvin mode reaches up to one fourth of the rotation velocity of the cylindrical container 
confirming that precession provides a rather efficient flow driving mechanism even at moderate 
values of Po. 

More relevant for dynamo action might be free Kelvin modes with higher azimuthal wave 
number. These free Kelvin modes are triggered by non-linear interactions and may constitute a 
triadic resonance with the fundamental forced mode when the height of the container matches 
their axial wave lengths. Our simulations reveal triadic resonances at aspect ratios close to those 
predicted by the linear theory except around the primary resonance of the forced mode. In that 
regime we still identify various free Kelvin modes, however, all of them exhibit a retrograde drift 
around the symmetry axis of the cylinder and none of them can be assigned to a triadic resonance. 
The amplitudes of the free Kelvin modes always remain below the forced mode but may reach up 
to 6% of the of the container’s angular velocity. The properties of the free Kelvin modes, namely 
their amplitude and their frequency, will be used in future simulations of the magnetic induction 
equation to investigate their ability to provide for dynamo action. 
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Instabilities and waves in rotating fluids are important in numerous technical 
applications. In most cases it is important to avoid these phenomena in order to 
ensure the stability of fast spinning liquid filled bodies like gyroscopes (Pfeiffer 1974), 
spacecrafts (Bao & Pascal 1997) or projectiles (Miller 1982). Waves in rotating fluids 
are also of general interest because of their fundamental character in oceanographic and 
atmospheric flows (Rossby 1939) like, e.g., zonal flows in the atmospheres of Saturn 
or Jupiter that may develop via a modulation instability of meridional Rossby waves 
(Connaughton et ah 2010). Moreover, the identification of inertial waves in the Earth’s 
fluid outer core (Aldridge & Lumb 1987) allows conclusions on the core dynamics and the 
associated dynamo process that is responsible for the generation of the Earth’s magnetic 
held. Nowadays, it is believed that in general planetary dynamos are driven by rotating 
thermal or compositional convection (King et ah 2010, M. Brener et ah 2010, D. Brener 
et ah 2010). On the other hand, other types of flows cannot be ruled out, and, for 
example, a superposition of random inertial waves in a rotating conducting fluid is 
capable of transferring energy in a magnetic field as well (Moffatt 1970). However, in 
the model of Moffatt (1970), held and flow must decay in the long term because of the 
lack of energy sources required for a steady driving. This problem can be overcome by 
more complex models in which, for example, the equatorial heat flux in the Earth’s outer 
core provides a persistent excitation mechanism for inertial waves with sufficient helicity 
for the generation of planetary magnetic fields (Davidson 2014). Inertial waves can be 
excited by precessional forcing as well, and it has long been discussed whether precession 
of the Earth can provide the necessary power to drive the geodynamo (Stewartson & 
Roberts 1963, Busse 1968, Malkus 1968). 

In a next-generation dynamo experiment currently under development at 
Helmholtz-Zentrum Dresden-Rossendorf (HZDR) a precession driven flow of liquid 
sodium will be used to validate its suitability to drive a dynamo (Leorat et ah 2003, 
Leorat 2006, Stefani et ah 2012, Stefani et al. 2015). Since the precessional forcing 
provides a natural driving mechanism, it represents a qualitatively new approach 
compared to previous dynamo experiments, where an optimized flow was driven either 
by impellers (Gailitis et al. 2000, Monchaux et al. 2007) or by a system of electromagnetic 
pumps (Stieglitz & Muller 2001). 

A small-scale water experiment is currently running at HZDR in order to examine 
the hydrodynamic properties of a precession driven flow in a cylindrical container with 
strong forcing and a large precession angle. Measurements of pressure fluctuations and of 
axial velocity profiles show three distinct flow regimes in dependence of the precession 
ratio: a laminar state (with only the forced mode) followed by a non-linear regime 
with various modes superimposed on the fundamental mode, eventually leading to a 
permanent chaotic state when a critical precession ratio is exceeded (Herault et al. 2015). 
Presumably, the ability of the flow to drive a dynamo will be different in the three 
regimes and although the precessional dynamo experiment at HZDR will allow magnetic 
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Reynolds numbers of the order of Rm ^ 700 (based on the rotation velocity and the 
radius of the cylindrical container) it is a priori not obvious whether a magnetic held 
will be self-generated at all. There are promising indications for dynamo action driven 
by precession both from liquid sodium experiments by Cans (1971) which - despite of 
the rather small precession ratio - achieved an amplihcation of an applied held by a 
factor of three, as well as from simulations in spheres (Tilgner 2005), spheroids (Wu 
& Roberts 2009), ellipsoids (Ernst-Hullermann et al. 2011), cubes (Krauze 2010) and 
cylinders (Nore et al. 2011). However, kinematic dynamo simulations demonstrate that 
the primary forced mode alone, which consists of a single, non-axisymmetric mode with 
azimuthal wave number m = 1 and axial wave number /c = 1, cannot drive a dynamo 
at a reasonable Rm (Giesecke et al. 2015). This coincides with Herreman & Lesahre 
(2011) who showed that a linear wave packet of inertial waves cannot drive a mean-held 
dynamo. However, this does not generally preclude the possibility of quasi-laminar (in 
the sense of Dudley & James 1989) or small scale dynamo action driven by inertial waves 
or another precession induced instability with sufficient intricate topology. 

Experimental studies in a weakly precessing cylinder hlled with a non-conducting 
huid showed a large diversity of how phenomena with complex structures frequently 
resulting in a vigorous turbulent state (Cans 1970, Manasseh 1992, Manasseh 1994, 
Kobine 1995, Manasseh 1996). Other experiments revealed large scale structures like 
a system of intermittent, cyclonic vortices that may provide a strong source of helicity 
and hence be benehcial for a dynamo (Mouhali et al. 2012), or free Kelvin modes with 
an azimuthal wave number m = 5 and m = 6 propagating around the symmetry axis of 
a weakly precessing cylinder (Lagrange et al. 2008, Lagrange et al. 2011). 

Free Kelvin modes are the natural eigenmodes in a rotating cylinder (Thomson 
1880, Greenspan 1968, McEwan 1970) and are promising candidates for driving quasi- 
laminar dynamo action in case of processional forcing because at least a subclass 
of these modes has a structure similar to the columnar convection cells that are 
responsible for dynamo action in convection driven models of the geodynamo (Lorenzani 
& Tilgner 2003). Free Kelvin modes may emerge from non-linear interactions of a forced 
Kelvin mode with itself or as a parametric instability that involves the forced mode and 
two free Kelvin modes (Kerswell 1999, Lagrange et al. 2008, Lagrange et al. 2011). 
The latter case requires appropriate combinations of wave numbers, frequencies and a 
geometry such that all three modes become resonant simultaneously which usually is 
named a triadic resonance. Triadic resonances have been observed in simulations of 
precession driven flow in a spheroid (Lorenzani & Tilgner 2003) and, experimentally, in 
weakly forced precession in a cylinder (Lagrange et al. 2008) and a cylindrical annulus 
(Lin et al. 2014). 

In the present study we conduct numerical simulations of a precessing flow far below 
the transition to the chaotic state observed in water experiments, but with sufficient 
forcing to excite free Kelvin modes with azimuthal wave numbers m > 1. We focus on 
the impact of the aspect ratio which - at least in the linear approximation - determines 
whether an inertial wave becomes resonant. We identify different free Kelvin modes by 
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Figure 1. Sketch of the problem setup. The cylinder rotates with flc and the rotation axis 
processes with ilp. In the simulations we varied the aspect ratio T = H/R and fixed the angle 
between rotation axis and precession axis at a = tt/2. 


means of their spatial structure, and we compute their frequencies from the azimuthal 
drift motion in order to conclude whether the free Kelvin modes constitute a triadic 
resonance. The results are used to further constrain the modes that may be observed 
in the water experiment at HZDR[|] and will be applied in future simulations of the 
magnetic induction equation in order to proof whether free Kelvin modes are suitable 
to drive a dynamo. 

2. Theoretical background 

The flow of a fluid with kinematic viscosity z/ in the frame of an enclosed precessing 
cylinder is described by the Navier-Stokes equation including terms for the Coriolis force 
and the Poincare force (see e.g. Tilgner 1998): 

du 

— + u ■ Vu + 2(f2p + flc) X u = —VP + z/V^u — (fip X flc) x r. (1) 

Here, u denotes the velocity held, which additionally obeys the incompressibility 
condition V ■ u = 0, r is the position vector, P is the reduced pressure (including 
the centrifugal terms), flc is the angular frequency of the cylindrical container and fip 
denotes the precession (see hgure [^. The precession axis is time dependent in the 
cylinder frame, and the unit vector fip that denotes the orientation of the precession 

f Pressure measurements in the water experiment at HZDR show a periodic signal with two frequencies 
close to the values expected from the dispersion relation for free Kelvin modes (see below). However, 
the measurements do not yet allow a unique identification of azimuthal, axial or radial wave numbers. 
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= sin a cos VLctx — sin a sin VLcty + cos az (2) 

with a the angle between rotation axis and precession axis. We consider a cylinder with 
aspect ratio F = H/R, where H and R are height and radius so that r G [0; i?] and 
z G [0; Fi?]. In the following we non-dimensionalize all quantities using R as length scale 
and 12“^ as time scale. 

In the inviscid limit a linear solution that fulhlls the boundary conditions z ■ u = Q 
at the top and the bottom is given by a Kelvin mode (Thomson 1880, Greenspan 1968): 


2:, (p, t) 




+ c.c. 


( 3 ) 


with a frequency ^k,m,n that depends on the azimuthal wave number m, the axial wave 
number k and a third number n that counts the roots of the dispersion relation 


GjAj Jm-i(Aj) + m (2 — Gj) Jmi^j) = 0 with Gj 



( 4 ) 


in which Jm denotes the Bessel function of order m and the triple index (fc, m, n) is 
replaced by j. Positive (negative) frequencies correspond to retrograde (prograde) 
propagation. Both signs correspond to different radial wave numbers (different solutions 
of the dispersion relation) and hence to a different spatial structure of the corresponding 
Kelvin mode. 

The dispersion relation ensures the fulhllment of the radial boundary conditions 
r ■ w = 0 at the sidewalls and although \j is not an integer it plays a role similar to 
a radial wave number with its position in the sequence of zeros corresponding to the 
number of half-cycles in the radial direction. The spatial structure of a Kelvin mode in 
a meridional plane Uj (r, z) is given by 


u^Jr,z) = 


u'j (r, z) = 


—i 


(4 - m 


Gj Aj (AjT*) -j- 


m(2 — Gj) 




(4 - 

u‘(r,z) = - i(^J„(A,r)sm ( 


2AV,„_i(A,r) - j„(A,r) 


COS 


cos 


/ Tlkz 
/ Tlkz 


( 5 ) 


Tlkz 


with Xj and XIj taken from the solutions of the dispersion relation (|^. 

In the inviscid linear approximation the total flow of the forced mode with m = 1 
is given by 


Um=i{r, <f) = X] X] AiknUiknir, (6) 

k n 

where we replaced the index j by the individual contributions and Ai^n represents the 
amplitude of the forced mode which can be computed by a projection of the applied 
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forcing onto a particular mode (Manasseh 1994). For a precession ratio (or Poincare 
number) Po = flp/flc and a precession angle a the corresponding calculation yields 
(Liao & Zhang 2012) 


_ -Posinar3(l - (-1)^)(2 + fij)2(2 - 12^) 
" “ A2[(ta/r)2 +1(2 - fij)](i - 


(7) 


If the wave length that corresponds to an axial wave number of a mode matches 
exactly the height of the cylinder, this mode becomes resonant with an eigenfrequency 
Vtj = 1 and the expression for the inviscid amplitude (j^ diverges. For each F there exist, 
in principle, an infinite number of resonant modes. From the dispersion relation (|^ we 
find that the primary forced mode, i.e. the mode with m = 1, the axial wave number 
k = 1 and a radial wave number corresponding to the first root, becomes resonant 
at F = 1.9898174 which is rather close to the so called spherical cylinder with height 
equals diameter (F = 2). Further resonances for increasing radial wave number occur 
at F = 0.9560735 (for n = 2) and at F = 0.6206981 (for n = 3). The computation of 
the amplitude at resonance requires a consideration of viscous effects in the bulk and in 
terms of boundary layers with associated Ekman layer suction. Then the amplitude can 
be calculated by matching the Ekman layer suction to the precessional force (Cans 1970). 
A more general approach for the computation of the amplitudes of the response created 
by the precessional forcing is given by Liao & Zhang (2012) who derived an asymptotic 
solution including viscous effects that is valid at and away from resonance. 

The precessional forcing only excites modes with m = 1 and k odd. Higher 
azimuthal modes or modes with even axial wave number must be triggered by non¬ 
linear interactions, e.g., in terms of triadic resonances that involve this forced Kelvin 
mode and two free Kelvin modes (ma, ka) and {mb, kb) or by the non-linear interaction 
of the forced Kelvin mode with itself (Meunier et al. 2008). In the following we assume 
that the resonant case is most promising in view of the dynamo problem due to the 
expected larger amplitudes of the velocity field. 

The free Kelvin modes are solutions of the linearized Navier-Stokes equation with 
a Coriolis term as a restoring force but without a precessional driving on the right hand 
side: 

3u 

— -h 2(f2p -h fic) X u = -VP. (8) 


The principle of the interaction of free and forced modes becomes evident when we 
consider an exact triad with a forced mode U{, with mf = 1 and fcf = 1 and two free 
Kelvin modes Ua with and kg, and Ub with and fcb- We ignore viscosity and 
further assume a constant amplitude Af of the forced mode so that we can write the 
total flow as 


U = AiUi + a{t)Ua + b{t)Ub, 


(9) 


where a{t) and h{t) denote the amplitude modulations of the free mode (without the 
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harmonic oscillations ^2^ and which are given by Kerswell (1999) 

^ = {Ua, t/b X (V X u;) + u;x{vx u^)) b{t), 

db 

- = ([/,, [/, X (V X C/f + t/f X (V X Ua)) ait). 

The square brackets {•, •) denote a scalar product dehned as U}y) = / C/* ■ U^dV, 
and a normalization of the free modes is assumed such that {Ui,Uj) = 6ij. In order 
to achieve non-trivial solutions of (10), and to constitute a triadic resonance, the wave 
numbers and frequencies of the free modes a and b must fulhll the conditions 


6m = mb — ma = m^ 

6k = \kb — ka\ = h 

6n = Tib — = Tic 


1 , 

1 , 

1 , 


( 11 ) 


where mf = 1, fcf = 1 and lie = 1 denote the properties of the forced mode. In contrast 
to sine and cosine-functions, there are no corresponding addition theorems for Bessel 
functions so that no further restrictions are imposed on the interaction of different radial 
wave numbers. 

In the following, we solve the dispersion relation (|^ for fc = 1, 2 and m = 1,..., 7. 
We assume that higher radial modes are damped and only consider the hrst root for 
each combination of m and k. Figure shows the frequency difference of two free 
Kelvin modes with 6m = 1 and 6k = 1 versus the aspect ratio in the vicinity of the 
critical value = 1. The crossings of the curves with the horizontal line at = 1 
mark aspect ratios for which an exact triadic resonance may occur. These aspect ratios 
at which 512 = 1 are listed in table together with the corresponding frequencies of the 
resonant free Kelvin modes. 



aspect ratio (height/radius) 

Figure 2. Frequency difference 6ft in dependence of the aspect ratio F for two free Kelvin modes 
which fulfill 6m = 1 and 6k = 1. The crossings of the curves with the horizontal black line at 
6ft = 1 mark the aspect ratio at which triadic resonances can be expected. 
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rria 


rUb 

kb 


klb 

P 

6 

1 

7 

2 

0.3825 

-0.6175 

1.67420 

7 

1 

6 

2 

0.3356 

-0.6644 

1.70478 

7 

2 

6 

1 

0.6494 

-0.3506 

1.72568 

6 

2 

7 

1 

0.6956 

-0.3045 

1.80074 

5 

1 

6 

2 

0.3888 

-0.6113 

1.87742 

6 

1 

5 

2 

0.3352 

-0.6648 

1.91146 

6 

2 

5 

1 

0.6489 

-0.3511 

1.94082 

5 

2 

6 

1 

0.7013 

-0.2987 

2.04266 

4 

1 

5 

2 

0.3974 

-0.6026 

2.14270 

5 

1 

4 

2 

0.3350 

-0.6649 

2.17832 

5 

2 

4 

1 

0.6485 

-0.3515 

2.22216 

4 

2 

5 

1 

0.7092 

-0.2907 

2.36888 


Table 1. Possible triadic resonances with 5m = 1, Jfc = 1 and (5r2 = 1 and corresponding aspect 
ratios. Negative frequencies present pro-grade modes and positive frequencies denote retrograde 
modes. Only the first radial wave number has been considered for each value of m and k. 


A triad is always formed by free Kelvin modes consisting of one prograde and one 
retrograde mode (Waleffe 1993). Note that the frequencies of free Kelvin modes that 
constitute different triads often are close (if we restrict to triads with /Cq = 1 and k}, = 2, 
see table making it difficult to reliably identify individual Kelvin modes in realistic 
setups with potential frequency shifts due to viscous and/or non-linear effects. 

3. Numerical model 

We conduct numerical simulations in a cylindrical geometry with the code SEMTEX 
that applies a spectral element Fourier approach for the numerical solution of the Navier- 
Stokes equation (Blackburn & Sherwin 2004). In order to simplify the implementation 
of the precessional forcing we switch to the precessional frame in which the precession 
is stationary and the cylinder rotates with a frequency hlc- Thus only a term for the 
Coriolis force appears in the Navier-Stokes equation which now reads 

du 

— + u ■ Vu -b 2f2p X u = —VP -b 1'V‘^u. (12) 

Due to the rotation of the cylinder in the precessing frame the boundary conditions 
for the azimuthal flow change to = rVtc (at endcaps and sidewall) whereas no-slip 
conditions are still applied for and Ur- 

The problem is described by four parameters, the Reynolds number dehned with 
the angular velocity of the container Re = DcP^/z/, the precession ratio (or Poincare 
number) Po = Dp/Dc, the aspect ratio P = H/R and the precession angle a with 
cos a = fip ■ fic/1 rip 11 riel- In the present study, we keep the precession ratio fixed at 
Po = 0.0141 and the precession angle is set to a = 7r/2. All simulations are performed at 
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Re = 6500 except for one run at Re = 10000 to briefly examine the impact of increasing 
Re. The aspect ratio is varied in the range T G [1.7; 2.4] in order to hnd the container 
geometry for which individual modes or triads become resonant. 

4. Results 

4 . 1 . Pattern of the total flow 

The simulations are started from an initial state with pure solid-body rotation u = re^. 
After switching on the Coriolis force at f = 0 the fluid shows a direct response in form of 
a large scale mode with m = 1. In dependence of the aspect ratio this is only a transient 
state until after a hnite time (which may take up to 1000 rotation periods or more) a 
quasi-steady state is reached. The typical flow pattern for T = 2 and Re = 6500 is 
presented in hgure where the nested isosurfaces show a snapshot of the axial velocity 
Uz at 30% (60%, 90%) of its maximum value. The flow is clearly dominated by a 




Figure 3. Volume rendering of the axial velocity Uz at Re = 6500 and F = 
2. The nested isosurfaces show Uz at 30% (60%, 90%) of its maximum value. An 
animation of the flow’s behavior in the statistically stationary state is available at 
https://www.hzdr.de/db/VideoDl?p0id=45098. 


velocity mode with m = 1 and k = 1. This is the primary Kelvin mode driven by the 
precessional forcing. Superimposed on the m = 1 mode we see some non-axisymmetric 
(time-dependent) disturbances that will be analyzed below. 

4-2. Kinetic energy 

4-2.1. Forced mode In the following we discuss the behavior of the flow by means of 
azimuthal Fourier modes Ujn{r,z) with the total flow given by u{r,z,ip) = -f 
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1) Y.m=o ^rn{r, . Figure 0 shows the history of the kinetic energy of the m = 1 
mode for different aspect ratios. A hnal state arises either after a short transitional 
period (F = 2.000) or after one or more bifurcations (see e.g. the blue curve for 
F = 2.150). The hnal state can be chaotic (F = 2.000, black curve), oscillatory 
(F = 2.15, blue), quasi-periodic with collapses (F = 2.2, red), quasi-periodic with bursts 
(F = 2.225, green) or stationary (F = 1.8, yellow). We also hnd transient periods with 
constant energy (e.g. F = 2.15 for t = 300...600), and it may take up to 1200 revolutions 
of the cylinder (F = 2.15, blue curve) till a hnal quasi-stationary state is reached. This 
saturation time also depends on the Reynolds number and decreases with increasing Re. 

The bifurcations shown by the m = 1 mode go along with the emergence of non- 
axisymmetric modes with m > 2 (hgurej^. These modes become unstable only after 
quite a long time, e.g., for F = 1.825 it takes up to 600 rotation periods of the cylinder 
till modes with m > 2 start to grow, and it may take up to twice that time till a quasi¬ 
steady state is reached. The non-axisymmetric modes with m > 2 essentially proceed 
similar to the behavior of the m = 1 mode, i.e., we see a chaotic behavior when the 
m = 1 mode behaves chaotic (F = 2.0, hgure central panel) and we see periodic 
growth and decay when the m = 1 mode is periodic (F = 1.825, top panel in hgure 
and F = 2.200, bottom panel in hgure |^. In the periodic regime the amplitude of the 
energy oscillations of the modes with m > 2 is rather large so that the kinetic energy 
may vary in time by more than one order of magnitude. 

In the following we compare analytic expressions for the energy of the forced mode 
with the time-averaged energy taken from the quasi-stationary period of the hnal state 
in our simulations. This comparison is only partly justihed because the analytical 
expressions are based on a number of serious simplihcations, such as the linearization of 
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Figure 4. Kinetic energy of the m = 1 mode versus time for various F. Note the qualitative 
diversity for different aspect ratios. 
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Figure 5. Temporal behavior of the modes with m = From top to bottom: 

r = 1.825,2.000,2.200. Note the different vertical scale for F = 2.000 (central panel). 


the Navier-Stokes equation, the limit Po —)■ 0, and the neglect of the time dependent part 
of the Coriolis force which is justihed only for small precession angles a. Furthermore, 
the hnal state in our simulations is instable and additional caution is advised when 
comparing simulations with a linear time independent solution. Nevertheless, such an 
analysis is helpful to identify the regimes in which, e.g., a maximum response of the 
flow can be expected, or to localize the aspect ratios at which triadic resonances are 
possible. 

Analytically, the energy of the forced Kelvin mode is given by (Liao & Zhang 2012) 


i?^=i = r2^ 


(Wr)^ + |(2-f^,) 




^ 1 ) 


|A,|Vi^(A,) 


(13) 


where the index j now denotes a combination of the axial wave number k and the radial 
wave number n, and Aj is the inviscid linear amplitude given by ([^ with the azimuthal 
wave number hxed at m = 1. Aj can also be calculated by the method of Liao & 
Zhang (2012) which includes viscosity and is valid at and off the resonance. However, in 
that case equation (13) is accurate only up to an error of 0 (Re“ 2 ) since equation (13) 


does not consider the energy from the flow in the boundary layers. Figure shows 
the time-averaged energy taken from the simulations (black curve) compared with the 
analytic solutions (green curve, inviscid solution and blue curve, viscous solution). For 


the calculation of the analytical solutions equation (13) is truncated at n = 6 and k = 6 


which is sufficient to reach convergence (off-resonance in the inviscid case) and avoids 
the occurrence of too many resonances in the inviscid case (green curve) which anyway 
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Figure 6. Time-averaged kinetic energy of the fundamental mode (m = 1) versus aspect ratio. 
Re = 6500, Po = 0.014, a = 7r/2. The dashed red lines denote the resonances with larger k and/or 
n. For the computation of the energy, equation (131 has been truncated at A: = 6 and n = 6 
(with prograde and retrograde modes respectively). Note that the inviscid solution (green curve) 
diverges at the resonances whereas the viscous solution has a finite maximum of « 0.021 at 
r « 1.9814. 


are suppressed in the more realistic viscous or non-linear computations. Comparing the 
numerical solutions (black curve) with the linear solutions we see signihcant deviations 
around the primary resonance (A; = l,n = 1, I/es ~ 1.9898) while outside of this 
regime the agreement of all curves remains good (say for T ^ 1.85 and for T ^ 2.05). 
Two key features characterize the behavior of the mean energy in our simulations. On 
the one hand we see a clear shift of the maximum to smaller aspect ratios and the 
maximum energy arises rather close to the aspect ratio expected for the resonance with 
k = 3 and n = 3. However, this resonance is suppressed by viscous effects and we 
do not see strong contributions with k = 3 and/or n = 3 in our simulations (see also 
section 4.5.1). Therefore the correlation is either coincidental or we have to assume an 
indirect impact of the {k = 3, n = 3) resonance. The second scenario is supported by 
the abrupt transition from the absolute maximum to a rather low energy state right 
below r 1.871. This jump is also connected to a modihcation of the character of the 
m = 1 mode which changes its behavior below T cs 1.86 from the chaotic regime to the 
periodic regime. Abrupt transitions of the energy of higher azimuthal modes can also be 
seen at the other (theoretical) resonances of the forced mode {k = 5,n = 5 and, a little 
further away /c = 5, n = 4, see hgure which will be revisited in the next paragraph. 
Beside the shift of the maximum to smaller aspect ratios the energy obtained in the 
simulations remains considerably smaller compared to the linear viscous solution. The 
maximum of the energy based on the amplitude of Liao & Zhang (2012) is 0.021 
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(at r 1.9814) which is roughly two times larger than the maximum kinetic energy 
obtained in our simulations ~ 0.0093 at F = 1.871). The differences are not 

surprising because the prerequisites for the computation of the amplitudes Aj are not 
well met in our simulations. In particular, we observe a significant impact of non¬ 
linear interactions in the simulations in terms of higher non-axisymmetric modes and 
the forming of an azimuthal shear flow that both draw energy from the forced m = 1 
mode. 

4-2.2. Higher azimuthal wave numbers We only find higher non-axisymmetric modes 
with a noteworthy amount of energy between F 1.81 and F 2.24. This corresponds 
approximately to the regime in which the energy of the m = 1 mode deviates from 
the linear prediction (see previous paragraph). For the sake of brevity we will call 
this regime the non-linear regime, assuming that the higher modes are triggered by 
non-linear effects. The kinetic energy of the modes with m > 2 is roughly 2 orders of 
magnitude smaller than the energy of the forced mode but essential characteristics such 
as the maximum at F = 1.871 and the abrupt drop of kinetic energy below this threshold 
are also reflected in the behavior of the higher non-axisymmetric modes (figure]^. The 



1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4 

aspect ratio 


Figure 7. Time-averaged kinetic energy of azimuthal modes with m > 2. Re = 6500, Po = 
0.014, a = 7r/2. The vertical lines indicate the aspect ratios for which triadic resonances could be 
expected according to the dispersion relation The labels denote the corresponding azimuthal 
and axial wave numbers with — (-I-) marking retrograde (prograde) modes. 


sudden drop of the energy for F < 1.871 pertains to all modes which indicates that this 
breakdown might be a global flow property that is not connected to an individual mode. 

A similar steep decrease of energy can also be found on the two outer edges of the 
non-linear regime. It is striking that all the abrupt transitions are close to theoretical 
resonances of the forced mode {k = 3,n = 3 ^ Fres ~ 1.8621, A; = 5,n = 5 ^ 
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Fres ~ 1.8142, and with some greater distance /c = 5,n = 4 Fres ~ 2.2911 , see 
figure]^. These resonances are, however, suppressed and cannot be directly observed 
in the simulations. So far, we cannot conclude whether the correlation of the energy 
jumps and the occurrence of (theoretical) resonances is a coincidence, or whether they 
may be some systematic relationship. 

Away from the primary resonance of the forced mode we find two further regimes 
with local maxima of the higher modes. For 1.80 ^ F ^ 1.85 we see a narrow window 
with the maximum around F 1.825 in which the m = 6 and m = 7 mode dominate 
(blue and yellow curve in figure]^. Likewise, but more explicit, we find the m = 4 and 
m = 5 mode becoming dominant for 2.05 ^ F ^ 2.25 with a maximum around F 2.2 
(green and red curves in figure . 

The local maxima around F 1.825 and F 2.2 for Fourier modes with 
5m = 1 indicate the presence of corresponding triadic resonances and roughly agree 
with predictions from the dispersion relation (see tableand figure]^. From the linear 
predictions we would also expect triads with m = 5 and m = 6 in the intermediate 
regime (for 1.85 < F < 2.05) but obviously there are other flow contributions with at 
least comparable energy which may disguise the energetic signature of the involved free 
Kelvin modes (e.g., we see relatively high energies of the m = 4 and the m = 5 modes, 
whereas in comparison the m = 6 mode is clearly suppressed). 

4 . 3 . Spatial Structure of the Fourier modes 

A qualitative impression of the flow structure is provided by isosurfaces of the individual 
Fourier modes of the axial velocity shown in figure for F = 2 (for an impression of 
the temporal behavior see the snapshots of the m = 5 mode in figure and the movie 
at https://www.hzdr.de/db/VideoDl?p0id=45104). The m = 1 mode is obviously 
dominated by a fc = 1 contribution and shows comparatively small variations in time. 
Regarding the higher non-axisymmetric modes we see a complex pattern for the m = 2 
and m = 3 mode with irregular temporal and spatial fluctuations and a more regular 
behavior of the modes m = 4, 5 and 6 with typical signatures of fc = 1 and k = 2. It 
is striking that the m = 6 mode, which should become resonant around F 2, remains 
weak compared to the m = 4 or the m = 5 mode. 

The amplitude of the individual contributions with m = 4, 5, 6 varies in time and 
the resulting flow perpetually changes its axial structure. The azimuthal mode resulting 
from the superposition of A; = 1 and k = 2 contributions has a distinct and time 
dependent asymmetry with respect to the equatorial plane of the cylinder, so that the 
associated flow is concentrated alternately in both halves of the cylindrical container 
(see time series of the m = 5 mode in figure]^. The typical timescale for this process 
is of the order of the rotation period 27r/r2c but we can not derive a reliable periodicity 
from our simulations. 

Away from the primary resonance of the forced mode the structure of higher 
azimuthal modes is much more regular and - except for periodic variations of amplitude 
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- essentially remains constant over time. Figure shows the dominant azimuthal 
modes (beyond m = 1) at F = 1.825 (m = 6 and m = 7) and at F = 2.200 (m = 4 and 



Figure 8. Snapshot of the axial velocity field Uz filtered to m = 1... 6 (upper left to lower 
right) for F = 2 (isosurface at 50% of the respective maximum value). An animation of the flow 
pattern is available at https://www.hzdr.de/db/VideoDl?p0id=45104. The temporal behavior 
of the modes m = 2 and m = 3 is complex and chaotic. The modes m = 4, 5 and 6 show a more 
regular behavior essentially resulting from the superimposition of free Kelvin modes with fc = 1 
and k = 2. Re = 6500, Po = 0.014, a = 7r/2. 



Figure 9. Time series of the to = 5 mode that visualizes the axial symmetry breaking. The 
isosurfaces show u^^^ z{,r,ip,z) at 50% of its maximum value at various characteristic time 
steps (form upper left to lower right: t=2617.5418, 2617.8489, 2618.0845, 2618.1958, 2618.8738, 
2618.9599 (time in units of the rotation time); note that the time gaps are not equidistant). 
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m = 4 m=5 


Figure 10. Snapshots of the structure of the axial velocity field for different Fourier modes. Left: 
r = 1.825, dominated by to = 6, A; = 1 and to = 7, A = 2. Right: F = 2.200, dominated 
by TO = 4, A = 2 and to = 5,A = 1. In both cases these modes represent the dominant 
contributions of the non-axisymmetric modes with to > 2 (see figure 13 below and animation 
at https://www.hzdr.de/db/VideoDl?pDid=45105). 


m = 5). In both cases, the axial behavior of the modes shows clear indications of axial 
wave numbers k = 1 and/or k = 2 and we will show below that these modes indeed 
fulhll all conditions of equation ( [II| ) and thus constitute triadic resonances. In contrast 
to the behavior at F = 2, the geometric structure of the dominant modes exhibits only 
minor temporal variations except strong oscillations of the amplitude (see animation at 
https:// WWW .hzdr.de/db/VideoDl?p0id=45105). 


4-4- Radial dependence 

In order to quantify the visual observations made in the previous section and to 
estimate amplitude and frequency of individual Kelvin modes we use a discrete sine- 
transformation (DST) applied to each Fourier mode 




N, 

(r) = 

< 7=1 


I sm 


qkTT 

iV. + l 


(14) 


with Nz the number of points in axial direction. Figure [TT] shows radial prohles for the 
dominant contributions to the axial velocity at F = 1.825 (a) and at F = 2.200 (b) where 
we only present the three leading modes for each case. The solid thick curves present 
the time average and the bright colors in the background represent the variation of the 
radial prohles in time. For comparison, the thick dashed curves show the analytical 
behavior oc Jmi^j^) predicted by the linear in-viscid approximation (|^ with Aj the 
hrst root of the dispersion relation (|^. The differences between the numerical data 
and the analytical curves are surprisingly small with deviations mainly in the boundary 
layers which can be explained by the different boundary conditions for the simulations 
(no slip) and the in-viscid solution (free slip). Except for (considerable) changes in the 
amplitude the pattern of the radial prohles exhibits only minor temporal changes and 
shows no signihcant contributions from higher radial modes (see color coded contour 


plots in hgure 11 that show the variations of the radial prohles in time). 
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k=1, time average) 
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k=1, n=1, analytic) 
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Figure 11. Radial dependence of the axial velocity decomposed into the dominant Fourier- and 
axial modes according to equation (14). (a) F = 1.825, from left to right: {m = l,k = 1), 
{m = 6 , k = 1) and (m = 7, k = 2), (b) T = 2.200, from left to right: (m = l,fc = 1), 
(m = 4, fc = 2) and (m = 5,k = 1). The solid thick curves show the temporal average and 
the weakly saturated colored curves show radial profiles at different time steps representing the 
temporal variations. The thick dashed curves denote the theoretical profiles from the linear 
inviscid solutions (oc Jm(Ajr), see equation §)■ The colored contour plots below each profile 
show the corresponding variation of the radial profiles in time. 


The behavior becomes more complex at T = 2 at which we see many modes with 
different m and k with nearly the same amplitnde and clear changes in the radial 
strnctnre. Fignre 12 shows characteristic radial profiles for the mode {m = 5, k = 1) 


and their temporal behavior as a typical example. The amplitnde exhibits considerable 
finctnations and, nnlike in the previons cases, we see additional changes in the radial 
strnctnre. When the amplitude is around a (local) maximum, the radial velocity profile 
is in accordance with a wavenumber that corresponds to the first root of the dispersion 
relation and the radial profiles can nicely be described by a function oc J 5 (Air) (see 
red curves in figure [T^). However, when the amplitude is weak we see clear signatures 
of higher radial wave numbers (see figure 12 b and d). Even though the higher radial 
modes remain weak, they might represent a channel for the transfer of energy from the 
(periodically occurring) maxima of a mode into smaller scales. 
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Figure 12. Snapshots of radial profiles ^(r) with m = 5 and /c = 1 at F = 2. (a) snapshots 
taken at local maxima of the (m = 5, /c = 1) mode dominated by n = 1 (marked by the the 
vertical red lines in panel (c)). The blue curve represents the analytical behavior oc J 5 {Xjr) with 
Xj the first root of the dispersion relation, (b), (d): snapshots of the radial profiles that match 
to a higher radial wave number. The red curves denote the numerical data and the blue curves 
denote a fit oc J^lXjr) with Xj corresponding to the second, third and fourth root of the dispersion 
relation Q. The corresponding time steps are marked by the green lines in panel c. (c) temporal 
variations of the radial profiles and of the amplitude \u\^z\ dX r — 0.67. 
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4-5.1. Height equal diameter (F = 2) Figure 13 shows the temporal behavior of \u 


If I 

m^z I 


at r = 0.67 for m = 1 to m = 7 (from top to bottom) and for k = 1 to k = A (blue, 
red, green, yellow curve). The peak values of the amplitudes are additionally listed in 
table (F = 2) and in table (F = 1.825 and F = 2.2) together with the frequencies 
obtained from the time derivative of the azimuthal phase of each mode. Note that 
we only list modes with a regular behavior of the frequencies, i.e., contributions which 
exhibit a steady and unique azimuthal drift that allows a conclusive computation of 
the time derivative of the azimuthal phase. The results quantitatively complement the 
observations made in the previous paragraph. For F = 2 (central column and table 
we see a dominant forced mode with m = 1 and k = 1 and the contributions with k > 1 
are negligible (smaller by a factor of 20). The higher azimuthal modes (m = 2... 7) 
reach approximately 10% of the amplitude of the m = 1 mode with k = 1 and k = 2 
slightly prevailing over k = 3 and k = 4. A striking property of the flow at F = 2 is the 
common orientation of the azimuthal drift motion of all modes (see table [^. We only 


Re 

m 

k 

max 

^num 

f^theo 

6500 

1 

1 

0.2673 

-1.0000 

— 1.0000 (forced) 


3 

1 

0.0412 

+0.2911 

+0.3342 (for n = 2) 


3 

2 

0.0266 

+0.3154 

+0.3333 (for n = 5) 


4 

1 

0.0536 

+0.0357 

no corresponding mode 


4 

2 

0.0525 

+0.3593 

+0.3647 (for n = A) 


5 

1 

0.0491 

+0.2278 

+0.2030 (for n = 3) 


5 

2 

0.0578 

+0.4692 

+0.5157 (for n = 2) 


6 

1 

0.0392 

+0.2006 

+0.1868 (for n = 3) 


6 

2 

0.0391 

+0.6239 

+0.6309 (for n = 1) 

10000 

1 

1 

0.2775 

-1.0000 

— 1.0000 (forced) 


3 

1 

0.0329 

+0.3090 

+0.3342 (for n = 2) 


3 

2 

0.0302 

+0.3013 

+0.3333 (for n = 5) 


4 

1 

0.0562 

+0.1687 

+0.1532 (for n = 5) 


4 

2 

0.0420 

+0.3511 

+0.3647 (for n = A) 


5 

1 

0.0570 

+0.2001 

+0.2030 (for n = 3) 


5 

2 

0.0574 

+0.5168 

+0.5157 (for n = 2) 


6 

1 

0.0417 

+0.3452 

+0.3204 (for n = 1) 


6 

2 

0.0437 

+0.6410 

+0.6309 (for n = 1) 


Table 2. Amplitudes and frequencies from the simulations at F = 2. Only modes with well- 
defined drift behavior are listed. The columns denote (from left to right): Reynolds number, 
azimuthal wave number, axial wave number, maximum of the amplitude, azimuthal drift frequency 
from the numerical data, and theoretical frequency from the dispersion relation Q for larger radial 
wave numbers n that fit best to the observed frequencies. 








Triadic resonances in a precessing cylinder 


20 


0.0001 

0-1000 















0.0001 

0-1000 




.-dtfa,,, 

. . , . . . 





ifwiii 



1 i Ir ' 1 
'\ 1 


time/rototion time 


2640 2660 2680 

time/rotation time 


2700 1560 1580 1600 1620 1640 1660 1680 

tlme/rototion time 


Figure 13. Temporal behavior of the amplitude of the flow ^ at r = 0.67 decomposed into 
Fourier- and axial modes with m = 1...7 (from top to bottom) and k = 1...4 (blue, red, green, 
yellow curve). From left to right: F = 1.825, F = 2.0 and F = 2.2. 
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find modes with a retrograde azimuthal drift so no combination is possible that fulhlls 
the requirements for a triadic resonance (in all cases (fhl 7 ^ 1, see table |^. The presence 
of distinct frequency signals that fit to higher radial modes (which however are hardly 


evident in the radial profiles for the (m = 5, A; = 1) mode, see figure 12) shows that 
there must be further contributions with chaotic behavior which are probably dominant 
and having no regular azimuthal drift. 

We also performed simulations at a slightly larger Reynolds number Re = 10000 
(at r = 2) and found only minor changes with respect to the run at Re = 6500 except 
for the modes (m = 4, /c = 1) and (m = 6 , fc = 1) which have lower frequencies at 
Re = 6500. However, so far our data is not sufficient to explain the impact of Re on 
the drift frequencies or to establish a conclusive scaling towards more realistic Re that 
will be reached in the experiment. Hence, we refrain from any further discussion of the 
properties of the flow at larger Re. 


4-5.2. Triadic resonances atV = 1.825 and T = 2.200 The behavior of the amplitude 
confirms the modified characteristics of the flow when the forced mode is off-resonance 
(T = 1.825, left column in figure 13, or T = 2.2, right column in figure 13; see also 
table 1^. For T = 2.200 we find a well-defined single triadic resonance with unique 
wave numbers (m = 4, fc = 2 and m = 5, /c = 1) and frequencies {Qa = —0.6368 and 
Qb = 0.3632). These values fulfill the conditions for a triadic resonance and are quite 
close to the theoretical values obtained from the dispersion relation (|^ (see table [^. 
The emergence of a triadic resonance is less explicit at T = 1.825. However, at this 


aspect ratio we hnd even two triadic resonances (see figure 13) with (m, fc) = (6,1) 


and (7, 2) and a second resonance with (m, k) = (5,1) and ( 6 , 2) (with much weaker 
amplitude, see table |^. In all cases the amplitude of the free modes that constitute a 


T 

m 

k 

max \u^\ 

^num 

^theo 

1.825 

1 

1 

0.1113 

- 1.0000 

— 1.0000 (forced mode) 


6 

1 

0.0272 

-0.2979 

-0.3045 (at T = 1.80074) 


7 

2 

0.0382 

+0.7010 

+0.6956 (at T = 1.80074) 


5 

1 

0.0094 

-0.2596 

-0.3511 (at T = 1.94082) 


6 

2 

0.0144 

+0.7405 

+0.6489 (at T = 1.94082) 

2.200 

1 

1 

0.1095 

- 1.0000 

— 1.0000 (forced mode) 


4 

2 

0.0596 

-0.6368 

-0.6650 (at T = 2.17832) 


5 

1 

0.0423 

+0.3632 

+0.3350 (at T = 2.17832) 


Table 3. Amplitudes and frequencies from the simulations for the dominant modes at F = 1.825 
and r = 2.200. The columns denote (from left to right): aspect ratio, azimuthal wave number, 
axial wave number, maximum of the amplitude, azimuthal drift frequency from numerical data, 
and theoretical frequency from the dispersion relation 0 at the aspect ratio at which the linear 
inviscid approximation predicts a triadic resonance (always supposing a radial wave number 
corresponding to the first root of the dispersion relation which is confirmed in figure [TT|). 
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triadic resonance strongly oscillates with a maximum of up to 50% of the forced mode 
(r = 2.2, m = 4, /c = 2) which corresponds to roughly 6% of the angular velocity of the 
cylindrical container. 


4-6. Azimuthal shear flow 


The growth of the triads goes along with a growth of an axisymmetric azimuthal shear 


flow that is mostly geostrophic (see left panel in hgure 14). The axisymmetric mode 


emerges with a slight delay with respect to the free Kelvin modes which indicates 
a saturation process by a detuning of the resonance frequencies (see right panel in 


hgure 14). The induced axialsymmetric how component is negative on average and, 


hence, causes a breaking of the solid body rotation which increases with increasing 
precession ratio. In the extreme case the breaking of the induced axisymmetric how 
entirely cancels the solid body rotation, giving the impression that in the laboratory 
frame the cylindrical container is rotating around a standing huid (Herault et al. 2015). 
This ehect is closely connected to the abrupt transition to a chaotic how which takes 
place at a critical precession ratio (Herault et al. 2015, Kong et al. 2014). 



Figure 14. Left: Pattern of the induced azimuthal shear flow at P = 2.2. Right: Temporal 
development of the amplitude at r = 0.67 of the azimuthal shear flow (green curve), the free 
Kelvin modes (red and green curve) and the forced mode (yellow curve). 


5. Conclusions 

We performed numerical simulations of a precession driven how in cylindrical geometry. 
The simulations conhrm that the energy that can be injected via precessional forcing is 
very sensitive to the aspect ratio of the cylindrical container. Signihcant contributions 
of higher non-axisymmetric components appear only in a limited range of aspect ratios 
(for T G [1.81; 2.24]) and it turns out that in that regime the behavior of the forced mode 
(with m = 1) cannot be described by a linear theory. Nevertheless, in all cases the forced 
Kelvin mode (m = 1) dominates and free Kelvin modes with higher azimuthal wave 
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number appear as distortions of the forced mode. The maximum response of the flow 
takes place at Tmax ~ 1.871 which is remarkable far away from the resonance predicted 
by the linear in-viscid approximation (at T = 1.98982). At the resonance maximum our 
results yield an amplitude for the ; 2 -component of the flow of ~ 0.27, which ~ 

assuming isotropy - roughly agrees with the estimations of Leorat (2006). The amplitude 
of the forced mode decreases with increasing distance from the resonance maximum, and 
in the regimes with triadic resonances, we still observe values of ~ 0.10 (in terms 

of the azimuthal velocity of the container). In real units (assuming ■ lOHz 

and i? = 1 m) this would correspond to fluid velocities of 15 m/s (at resonance) and 
6 m/s (off resonance) which is of the same order as the typical flow speed in the Riga 
dynamo experiment. Regarding the higher Kelvin modes, the saturated amplitudes 
of the strongest free Kelvin modes are mostly independent of the forced mode and 
achieve values of about 0.05 to 0.06 in terms of the angular velocity of the container 
corresponding to ~ 3.5 m/s in the experiment. 

At the aspect ratio envisaged for the experiment (T = 2 with H = 2m. and i? = 1 m) 
we may expect an energy of the fluid flow at roughly two thirds of the maximum value 
at the optimum aspect ratio (T = 1.871). However, this optimum aspect ratio probably 
depends on the Reynolds number and on the forcing, so that we believe that the chosen 
geometry for the dynamo experiment is a good compromise to ensure an efficiently driven 
fluid flow in a container with a ratio of diameter to height that at least approximately 
reflects the geometry of planetary bodies. 

Although at present we do not know whether the flow helds that emerge in 
our simulations will be able to drive a dynamo, we may try to estimate typical 
magnetic held strengths that can be expected in the planned dynamo experiment. We 
assume a saturation of the magnetic energy at roughly 5% of the kinetic energy of 
the hydrodynamic how which is substantiated by the saturation behavior of the Riga 
dynamo (Gailitis et al. 2008). We further presuppose that the internal velocity caused 
by the precessional forcing linearly scales with the rotation of the container (and hence 
with Re) and assume an angular velocity of 0,^ = 2t[ ■ 10 Hz, a radius R = 1 m, a height 
H = 2m and a density of liquid sodium of p ~ 930 kg/m^ (at ~400 K). After the change 
over to real units the kinetic energy obtained in our simulations at T = 2 (Rkin ~ 0.006) 
corresponds approximately to a magnetic held strength B ^ A/R05po//EAA2f ~ 40 mT. 
This is in the range of the Riga dynamo which is not particularly surprising since the 
typical velocities obtained in our simulations indeed match the velocity arising in the 
Riga dynamo. However, given that the corresponding m = 1 mode has produced no 
dynamo in the kinematic simulations it might be more realistic to refer to the free Kelvin 
modes with higher m which have less kinetic energy. In that case the typical magnetic 
held strength amounts to only B ^ 5 mT which, however, is still in the range of the 
values obtained at the Von-Karman-Sodium dynamo (Monchaux et al. 2007). We expect 
a stronger response and hence a larger saturation held strength for increasing Po, at 
least as long as we remain below the critical Poincare number for the sudden transition 
to a chaotic state (Po“‘* ~ 0.0725 at Re = 5.65 x 10^, see Herault et al. 2015). However, 
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unless we know whether or not the precession driven flow will have the right structure 
and sufficient magnitude to excite dynamo action at all, these estimations are only useful 
to demonstrate that the reference values assumed for the planned dynamo experiment 
will allow magnetic fields with a reasonable field strength. 

Regarding the structure of the flow, we And various manifestations of non- 
axisymmetric contributions in terms of free Kelvin modes with higher azimuthal wave 
numbers m. We And free Kelvin modes in resonance with the forced mode only when the 
forced mode is sufficiently far away from its primary resonance, i.e., when the amplitude 
of the forced mode is not too strong. Triadic resonances can be clearly identified using 
a combined Fourier-Discrete Sine transformation in the azimuth and along the axis 
with structure and frequencies close to predictions from linear theory. Actually, we 
were also expecting triads around F 2 with m = 5 and m = 6. Indeed, we do find 
patterns of free Kelvin modes in that regime, but no combination of them satisfies the 
triadic resonance conditions. Instead, we see free Kelvin modes solely with retrograde 
drift motion and for each azimuthal wavenumber m the contributions with k = 1 and 
k = 2 are approximately of equal strength so that in sum we see a clear breaking of 
the equatorial symmetry of the flow. This symmetry breaking has been essential for the 
functioning of a dynamo in simulations of Tilgner (2005) in a precessing sphere whereas it 
seemed less important in the direct numerical simulations of Nore et al. (2011). However, 
preliminary kinematic dynamo simulations we conducted recently with an analytic flow 
of Kelvin modes suggest that the joint appearance of contributions with k = 1 and 
k = 2 considerably facilitates the occurrence of dynamo action. 

The triadic instability requires a rather long time to emerge and saturates to 
a final state with slow periodic growth and decay. This period is likely related to 
the precession time scale but a conclusive correlation requires further simulations 
at different Po. The periodic growth followed by a fast decay of the free Kelvin 
modes in the regime with triadic resonances reminds of the so called resonant collapse 
(McEwan 1970, Manasseh 1992). However, in our simulations the energy is essentially 
distributed among very few, large scale modes. These modes are the fundamental 
forced mode and two free Kelvin modes that establish the triad and an axisymmetric, 
essentially geostrophic mode. We do not spot any turbulence-like behavior during the 
decay of the free Kelvin modes, i.e. there is no transition from the dominance of the large 
scale mode into small scale turbulent flow which would be characteristic for the collapse 
phenomenon. This behavior might change when the Reynolds number approaches larger 
(i.e. more realistic) values. This is substantiated by the simulations at F = 2. Here we 
see no triadic resonance, but rather a quasi-periodic occurrence of free Kelvin modes 
that exhibit a trend towards larger radial wave numbers during their decay which may 
represent (part of) a cascade to small scale structures. 

In principle the occurrence of very few dominant modes calls for a low-dimensional 
model. However, the ansatz (|^ used in this study to explain the occurrence of triadic 
resonances is too simplistic to allow for a reasonable description of the non-linear 
behavior and the saturation obtained in the simulations presented above. A more 
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sophisticated approach is presented by Lagrange et al. (2011) who use a low-dimensional 
model with four coupled ordinary differential equations that describe the amplitude of 
the forced mode, the free Kelvin modes and an axisymmetric geostrophic mode. This 
model allows for saturation via the axisymmetric mode and for a reinforcement of the 
forced mode by the free Kelvin modes and qualitatively reproduces the results obtained 
in our simulations. However, in order to yield a quantitative agreement, advanced 
models are still required. 
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